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Abstract —This paper focuses on a class of important two-hop 
relay mobile ad hoc networks (MANETs) with limited-buffer 
constraint and any mobility model that leads to the uniform 
distribution of the locations of nodes in steady state, and develops 
a general theoretical framework for the end-to-end (E2E) delay 
modeling there. We first combine the theories of Eixed-Point, 
Quasi-Birth-and-Death process and embedded Markov chain to 
model the limiting distribution of the occupancy states of a relay 
buffer, and then apply the absorbing Markov chain theory to 
characterize the packet delivery process, such that a complete 
theoretical framework is developed for the E2E delay analysis. 
With the help of this framework, we derive a general and exact 
expression for the E2E delay based on the modeling of both 
packet queuing delay and delivery delay. To demonstrate the 
application of our framework, case studies are further provided 
under two network scenarios with different MAC protocols to 
show how the E2E delay can be analytically determined for a 
given network scenario. Einaliy, we present extensive simulation 
and numerical results to illustrate the efficiency of our delay 
analysis as well as the impacts of network parameters on delay 
performance. 

Index Terms —mobile ad hoc networks (MANETs), limited 
buffer, end-to-end delay, performance modeling 


I. Introduction 

With the development of wireless communication tech¬ 
nologies, mobile ad hoc networks (MANETs) have become 
an appealing candidate for many critical applications, such 
as emergency rescue, disaster relief, coverage extension for 
cellular networks, etc. m-a. Although lots of work has 
been done to facilitate the commercialization of MANETs, 
understanding their fundamental delay performance has been a 
critical research issue for them to support various applications 
with different quality of service (QoS) requirements 16], a. 

End-to-end (E2E) delay, the time that a packet takes to reach 
its destination after it is generated by its source, serves as 
the most fundamental delay metric. The available theoretical 
studies on E2E delay of MANETs mainly focus on deriving 
its upper bound or approximation. Regarding the delay upper 
bound of MANETs, Neely et al. lH derived some useful 
results for a cell-partitioned MANET with the two-hop relay 
(2HR) routing scheme and i.i.d mobility model. Later, Gamal 
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et al. M and Sharma et al. ifTOl extended the results of lO 
to the continuous network model and general mobility model, 
respectively. Inspired by these works, extensive research ac¬ 
tivities have been devoted to the study of delay upper bound 
for MANETs under various network scenarios, such as under 
the motioncast in ifTTIl . under the cognitive networks in ifT^ . 
under the packet redundancy in m, under the multi-hop 
back-pressure routing in m, and under the power control 
in M- Regarding the delay approximation, Jindal et al. ca 
explored recently the E2E delay approximation for MANETs 
with multi-hop relay routing, and Liu et al. studied the E2E 
delay approximation for MANETs with probing-based 2HR 
routing iflTl and limited packet redundancy ifTSll . 

In addition to the studies on delay upper bound or approxi¬ 
mation for MANETs, Neely et al. |8l also applied the queueing 
theory to derive the exact expression for E2E delay. Eollowing 
this line, recently some results have been reported on the 
modeling of really achievable E2E delay in MANETs. Chen 
et al. 1191 explored the MANETs with Aloha MAC protocol 
and determined the corresponding exact E2E delay there under 
the continuous network model. Eor a cell-partitioned MANET 
with broadcast-based routing scheme, Gao et al. proposed 
a new theoretical framework for the analysis of its exact E2E 
delay based on the theory of Quasi-Birth-and-Death process. 

It is notable that the common limitation of above studies 
is that to simplify their analysis of E2E delay, they assume 
the relay buffer of a node, which is used for temporarily 
storing packets of other nodes, has an infinite buffer size. 
In a practical MANET, however, the buffer size of a mobile 
node is usually limited due to both its storage space limitation 
and computing capability limitation. Thus, for the practical 
delay performance study of MANETs, the constraint on buffer 
space should be carefully addressed. Notice that the E2E 
delay modeling with practical limited-buffer constraint still 
remains a technical challenge. This is mainly due to the lack 
of a general theoretical framework to efficiently characterize 
the highly dynamic behaviors in such networks, like the 
complicated buffer occupancy states of a relay buffer, as well 
as the highly dynamic queuing process and delivery process 
of a packet. 

As a step towards the modeling of real achievable E2E delay 
for the practical MANETs with buffer constraint, we focus on 
a class of important 2HR MANETs with limited shared relay 
buffer and propose a general theoretical framework for the E2E 
delay modeling there. The main contributions of this paper are 
summarized as follows. 

• Eor the concerned MANET, we first combine the theo- 
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ries of Fixed-Point (FP), Quasi-Birth-and-Death (QBD) 
process and embedded Markov chain (EMC) to construct 
an analytical model to fully depict the complicated oc¬ 
cupancy behaviors of a relay buffer with limited buffer 
size. 

• Based on the above modeling of relay buffer occupancy 
behaviors, we then apply the absorbing Markov chain 
(AMC) theory to characterize the packet delivery process, 
such that a complete theoretical framework is developed 
for the E2E delay modeling in the concerned buffer- 
limited MANETs. This framework is general in the sense 
that it can be applied to conduct E2E delay analysis for 
a 2HR MANET with any mobility model that leads to 
the uniform distribution of the locations of nodes, such 
as the i.i.d mobility model ||8], the random walk model 
||9l , the random way-point model ll2Tll . etc.. 

• To demonstrate the application of the proposed frame¬ 
work, case studies are further provided under two net¬ 
work scenarios, i.e., the cell partitioned networks with 
local scheduling-based MAC protocol (LS-MAC) ||8l 
and Equivalent-Class based MAC protocol (EC-MAC) 
II 22 I . to show how the E2E delay can be analytically 
determined for a given network scenario by applying our 
framework. Einally, extensive simulation and numerical 
results are provided to validate the efficiency of the 
proposed E2E delay model and also to illustrate the 
impacts of network parameters on delay performance. 

The remainder of this paper is organized as follows. Sec¬ 
tion |II] introduces preliminaries involved in this paper. The 
complicated relay buffer occupancy behaviors are analyzed in 
Section |III] We derive the queuing delay, delivery delay and 
E2E delay in Section |IV] and conduct case studies in Sec¬ 
tion The simulation results and corresponding discussions 
are provided in Section |VI] Einally, Section Ivnl concludes this 
paper. 

II. Preliminaries 

In this section, we first present some basic assumptions and 
the buffer constraint, and then introduce the routing scheme 
and some critical definitions involved in this study. 

A. Basic Assumptions 

We consider the following minimal set of assumptions; 

(A.i) The ad hoc network is time-slotted and consists of n 
mobile nodes. 

(A.ii) The packet generating process in each source node is 
independent and assumed to be a Bernoulli process, 
where a packet is generated by its source node with 
probability A in a time slot. 

(A.iii) The widely-used permutation traffic model ||8l, ll23l . 
Il24l is adopted. With this traffic model, there are n 
unicast traffic flows in the network, each node is the 
source of one traffic flow and also the destination of 
another traffic flow. We denote by tf{i) the destination 
node of the traffic flow originated from node i, then 
the source-destination pairs are matched in a way 



that the sequence , p{2), ■ ■ ■ ,i^(n)} is just a 

derangement of the set of nodes {1, 2, • • • ,n}. 

(A.iv) During a time slot the total amount of data that can 
be transmitted from a transmitter to its corresponding 
receiver is fixed and normalized to one packet. 

(A.v) We consider the mobility model that leads to the 
uniform distribution of the locations of nodes in steady 
state, which covers many typical mobility models such 
as the i.i.d mobility model, the random walk model, 
the random way-point model, etc.. More formally, we 
denote by Xi{t) the location of ith node at time slot 
t and assume the process {2fi(-)} is stationary and 
ergodic with stationary distribution uniform on the 
network area; moreover, the trajectories of different 
nodes are independent and identically distributed. 

B. Buffer Constraint 

As illustrated in Eig. [T] each node in the MANET maintains 
n — 1 individual queues, one source queue for storing the 
packets that are locally generated at this node, and n — 2 
parallel relay queues for storing packets of other flows (one 
queue per flow). All these queues follow the EIEO (flrst-in- 
flrst-out) discipline. 

Similar to the available studies on buffer-limited wireless 
networks Ea, m, we consider the following practical buffer 
constraint that all the n — 2 relay queues of a node share 
a common relay buffer with the limited buffer size of B 
packets, while the buffer size of source queue is unlimited. We 
adopt this buffer constraint here mainly due to the following 
reasons. Eirst, the mathematical tractability of this assumption 
allows us to gain important insights into the structure of E2E 
delay analysis. Second, the analysis under this assumption 
provides a meaningful theoretical result in the limit of infinite 
source buffer. Third, in a practical wireless network, each node 
usually prefers to reserve a much larger buffer space for storing 
its own packets than that for storing packets of other flows. 
Also, even in the case that the buffer space of source queue is 
not enough when bursty traffic comes, the congestion control 
in the upper layer can be executed to avoid the loss of locally 
generated packets ll25l . 

C. Handshake-based 2HR Scheme 

Regarding the routing scheme, we focus on the 2HR 
scheme, because it is simple yet efficient and thus serves as a 
class of attractive routing protocols for MANETs IH], ||23l. To 
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avoid unnecessary packet loss and support the efficient opera¬ 
tion of the concerned buffer-limited MANETs, we introduce a 
handshake mechanism with negligible overheacQ into the 2HR 
scheme such that the packet dropping will not happen even in 
the case of relay buffer overflow. Once a node (say S) gets 
access to the wireless channel in a time slot, it executes the 
new handshake-based 2HR (H2HR for short) routing scheme 
summarized in Algorithm [1] 


Algorithm 1 H2HR algorithm 
1: if The destination D is within the transmission range of 

S then 

2 : S executes Procedure [U 

3: else if There exist other nodes within the transmission 
range of S then 

4: With equal probability, S selects one node as the 

receiver. 

5: S executes Procedure |2] or Procedure [3] equally with 

the receiver. 

6 : end if 


Procedure 1 Source-to-destination (S-D) transmission 
1: if S has packets in its source queue then 
2: S transmits the head-of-line (HoL) packet in its source 

queue to D. 

3: S removes the HoL packet from its source queue. 

4: S moves ahead the remaining packets in its source 

queue. 

5: else 

6: S remains idle. 

7: end if 


Procedure 2 Source-to-relay (S-R) transmission 
1: if S has packets in its source queue then 
2 : S initiates a handshake with the receiver to check 

whether the relay buffer of receiver is full or not. 

3: if The relay buffer of receiver does not overflow then 

4: The receiver dynamically allocates a new buffer space 

to the end of the corresponding relay queue. 

5: S transmits the HoL packet in its source queue to the 

receiver. 

6: S removes the HoL packet from its source queue. 

7: S moves ahead the remaining packets in its source 

queue. 

8 : end if 

9: else 

10: S remains idle. 

11 : end if 


*The handshake mechanism can be easily implemented by sending only 
one indicator bit from the receiver to the transmitter (e.g., bit 0 when the 
relay buffer is full, and bit 1 otherwise), so the impact of this overhead can 
be neglected in our analysis. 


Procedure 3 Relay-to-destination (R-D) transmission 

1: if S has packets destined to the receiver then 

2: S transmits the HoL packet in its corvesponding relay 

queue to the receiver. 

3: S removes the HoL packet from this relay queue. 

4: S moves ahead the remaining packets in this relay 

queue. 

5: This relay queue releases one buffer space to the 

common relay buffer of S. 

6 : else 

7: S remains idle. 

8 : end if 


D. Definitions 

Here we introduce some important definitions involved in 
this study. 

Relay-buffer Overflowing Probability (ROP); Lor the 

concerned MANET with a given packet generating rate A in 
each node, the relay-buffer overflowing probability Po(A) of a 
node is defined as the probability that the relay buffer of this 
node overflows (i.e, the relay buffer is full). 

Queuing Delay: The queuing delay is defined as the time 
it takes a packet to move to HoL in the source queue (i.e., 
the source node starts to deliver it) after it is generated by its 
source. 

Delivery Delay: The delivery delay is defined as the time 
it takes a packet to reach its destination after its source starts 
to deliver it. 

End-to-end Delay: The end-to-end delay is defined as 
the time it takes a packet to reach its destination after it is 
generated by its source, which is the sum of its queuing delay 
and delivery delay. 

III. Relay Buffer Analysis 

In this section, we first introduce three basic probabilities. 
Based on these probabilities, we then apply the QBD process 
modeling and EMC technique to depict the occupancy behav¬ 
iors of a relay buffer. Linally, we construct a self-mapping 
function for the ROP Po(A) (i.e., Po(A) is the fixed-point of 
this function) to determine the limiting distribution of the 
occupancy states, which will help us conduct delay analysis 
in Section HVl 

Due to the symmetry of nodes and traffic flows, we only 
focus on one node S in the following analysis. We denote by 
Psd, Psr and Prd the probabilities that in a time slot S gets 
access to the wireless channel and decides to execute S-D, 
S-R and R-D transmission respectiveljH These probabilities 
can be determined under a given network scenario and the 
derivation of them will be elaborated in case studies. 

A. QBD Process Modeling 

Regarding the source queue of S, it can be modeled as a 
Bernoulli/Bernoulli queue lIZTll with packet arrival rate A and 

^It is notable that psr = Prd< ^nd executing a transmission doesn’t mean 
that S will successfully transmit a packet in this time slot. 
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Fig. 2. Transition cases from a general state 



Fig. 3. State transition diagram of the QBD process. 


service rate /is(A), where /is (A) is given by 

/is(A) =Psd+Psr(l-Po(A)). (1) 

Due to the reversibility of Bernoulli/Bernoulli queue, the 
packet departure process of source queue is also a Bernoulli 
process with rate A. 

Regarding the relay buffer of S, we adopt a two-tuple 
X(f) = to define its state at time slot t, where 

I {t) denotes the number of packets occupying the relay buffer, 
and J{t) denotes the number of relay queues which are not 
empty, here 0 < I{t) < B, I < J{t) < I{t) when I{t) > 0, 
and J{t) = 0 when I{t) = 0. 

As illustrated in Fig. |2l suppose that the relay buffer of S is 
in state {i,j) at the current time slot, only one of the following 
transitions may happen in the next time slot; 

Case 1: i < B,sl packet enters the relay buffer, and this packet 
is destined for a destination same as one of packet(s) 
already in relay queues. 

Case 2: i < B, a packet enters the relay buffer, and the des¬ 
tination of this packet is different from all packet(s) 
already in relay queues. 

Case 3: * > 0, a packet in one of the relay queues is delivered 
to its destination, and there still exist other packet(s) 
in this relay queue. 

Case 4: i > 0, a packet in one of the relay queues is delivered 
to its destination, and there is no remaining packet in 
this relay. 

Case 5; no packet enters into or departs from the relay queues. 

To facilitate our discussion, we call the subset of states 
Li = {(i,l),(i,2),--- level i, Lq = {(0,0)} level 

0, and state {i,j) that the relay buffer is in level i and 
phase j. Notice that when the relay buffer is in some state 
of level z at a time slot, the next state after one-step state 
transitions could only be some state in the same level or 
its adjacent levels. Thus, as time evolves, the state transi¬ 
tions of the relay buffer of S form a two-dimensional QBD 
process {X(f), f = 0,1, 2, • • • } Il28l . 1291 . According to the 
transition cases in Fig. |2] the overall transition diagram of 
the QBD process is summarized in Fig. [3] There are in total 
l-|-0.5i3(l-|-i3) states for the QBD process, and we arrange all 
these states in a low-to-high level and low-to-high phase way 



Fig. 4. State machine of the EMC. 


as follows: {(0, 0), (1,1), (2,1), (2, 2), • • • , (B, B)}. Then the 
corresponding state transition matrix P of the QBD process 
can be determined as 


Ao,o Aq i 

Ao,i Al l Ai_2 


Ab-1,B-2 Ab-1,B-1 Ab-i,b 
Ab,b-i Ab,b 


( 2 ) 


where the sub-matrix Ai.i is of size i x I (Aq.o, Aq i and 
Ai,o are of size 1 x 1), denoting the transition probabilities 
from the states of level i to the states of level 1. 

It is notable that in our QBD process of relay buffer, 
different levels have different number of phases, and the 
transition probabilities of one state depend on its level, thus the 
QBD process is level-dependent and it is very difficult to solve 
its limiting distribution by determining its critical matrices 
and conducting recursive algorithm 1281 . ll2^ . To address this 
issue, we adopt a Markov chain-collapsing technique lIMIl . 
ED to convert the two-dimensional QBD process to a one¬ 
dimensional EMC in the next subsection. 


B. Collapsing to an EMC 

For the QBD process of Fig. [2 we integrate all states of 
a level into only one state, then the two-dimensional QBD 
process is collapsed to a one-dimensional Embedded Markov 
Chain (EMC). As illustrated in Eig. |4l one state Li of the 
EMC corresponds to one level i of the QBD process, and 
Pj^’ denotes the one-step transition probability from state Li 
to state Li in the EMC. According to the EMC theory EOll . 
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EH, the state transition probability of the EMC is the phase- 
averaged state transition probability of the QBD process, then 
we have 


P(0,0),Lp 


i = 0 


{i,i) 

Pl = 


Pj\Li, ^<i<B 
i=i 


(3) 


where P(ij)^Li denotes the transition probability from state 
{i,j) to the states of level I, and Pj\Li denotes the conditional 
probability that the relay buffer is in phase j given that it is 
in level i. Based on formula Q as well as the ergodic and 
uniform features of the distribution of node location, we have 
the following lemma regarding the transition probabilities of 
the EMC. 

Lemma 1: The one-step transition probability of the 
EMC is determined as 


Ps ('^) ■ Psr 

i 


= { n-i + i 


’ Prdi 


1-Pl -Pl \ 


vO, 


I = i + 1 < B 

l=i-l>Q 

I = i 
others 


(4) 


where p^{\) = =---^ 

PsW Psd + Psr{l - Po{a)) 

Proof: See Appendix for the proof. ■ 

We arrange all the states of EMC in a low-to-high level 
way as follows: {Loi^ir’’ Then the corresponding 

state transition matrix P emc can be determined as 


P EMC 


( 0 , 0 ) 

pI 

(1.0) 

pI 


( 0 , 1 ) 

pI 

( 1 . 1 ) 

pI 


Pl 


1 , 2 ) 


(5) 


(B.B-l) (B.B) 

pI pI 


C. Constructing the self-mapping function 

With the help of transition matrix Pemc, we then construct 
a self-mapping function forpo(A), i.e., Po(A) is the fixed-point 
of this function ll^ . such that Po(A) as well as the limiting 
distribution of the occupancy states of a relay buffer can be 
determined. 

Erom the the state machine of EMC in Eig. |4] and the 
transition matrix Pemc^ we can see that: 1) the EMC is 
irreducible; 2) each state Li is recurrent; 3) the period of 
each state Li is 1, so each state is aperiodic. Based on these 
properties, we can conclude that the EMC is ergodic, thus its 
limiting distribution TEl = [Tr^g, • • • ^elb] exists and is 
unique, and is same as its stationary distribution ||33l. Then 
we have 


Hl • Pemc = IIl, (6) 

Hl • 1 = 1, (7) 

where 1 is a column vector of size (i7-|-l) x 1 with all elements 
being 1, and equation d?) follows from the normalization 


property of a probability vector. Combining (|6]l with (Q we 
have 


'^Li = 


c,-ps{\y 

EtoCk-Ps{xr' 


( 8 ) 


where Ci = 

It is notable that the relay buffer overflows when it is in 
level B, then the critical self-mapping function for Po(A) is 
constructed as 



PoW = f{Po{X)) = TTLb 


Cb-Ps{X)^ 

■ PsiX)^ 


(9) 


Given a packet generating rate A, the self-mapping function 
doesn’t contain any unknown parameters except Po(A). Thus 
by solving equation (|9ll, we can determine the ROP Po(A) 
corresponding to a given A, and the limiting distribution of 
the EMC can be recursively determined as 

'ttLi =Po{x) ■ ps{xy~^ ■ ( 10 ) 

'^B 

The limiting distribution 11 = [ttop, tti i, • • • , tt^ • • • , tte^b] 
of the QBD process can be further determined as 


= TtLi ■ Pj\Li, ( 11 ) 

where Pj\Li is given by formula (l42l) in Appendix lAl 

Remark 1: Notice that if we don’t apply the handshake 
mechanism, we can also develop the corresponding theoret¬ 
ical framework in the same way to model the relay buffer 
occupancy process, where the ROP Po{X) derived in (|9ll just 
corresponds to the packet dropping probability. 


IV. Delay Analysis 

With the help of ROP and limiting distribution of occupancy 
states of a relay buffer, in this section we analyze the delay 
performance for the concerned buffer-limited MANET. We 
denote by Q, D and T the queuing delay, delivery delay and 
E2E delay of a packet respectively. The E2E delay of a packet 
will be derived by computing its queuing delay and delivery 
delay respectively. The queuing delay will be obtained by 
analyzing the queuing process of the source queue, while the 
delivery delay will be derived by modeling the packet delivery 
process as an AMC and analyzing the time the chain takes to 
enter the absorbing state. 

Before presenting our main results on the delay perfor¬ 
mance, we first provide the following lemma regarding the 
per node throughput capacity, which is the maximal packet 
generating rate the MANET can stably support, and the 
corresponding delay can then be determined. 

Lemma 2: Eor the considered MANET, its per node 
throughput capacity p is given by 


P-Psd+Psr „ , „• 

n — 2 -\- B 

(12) 

Proof: See Appendix iBl for the proof. 

■ 
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A. Queuing Delay 

Considering a given packet generating rate A (A < /i), the 
corresponding ROP Po(A) can be obtained by solving equation 
and the service rate of source queue /is(A) can be further 
determined by formula ([TJ. Thus, in the following analysis, 
we use po and ps to represent Po(A) and /is(A) respectively 
if there is no ambiguous. Notice that the source queue is 
a Bernoulli/Bernoulli queue, thus its average queue length 
Lsource is given by ||22| 

— A — A^ 

^source — T- (13) 

Ms - A 

According to the Little’s Law ll34l . the average delay of a 
packet in its source queue E{7As} is given by 

nos} = (14) 

Ms - A 

Then, the expected queuing delay E{(5} is determined as 

E{Q} = E{iAj- —= . (15) 

Ms Ms(Ms-A) 


B. Delivery Delay and End-to-end Delay 

We present the following theorem regarding the expected 
E2E delay of the concerned buffer-limited MANET. 

Theorem 1: (Main result) Eor the concerned MANET with 
number of nodes n, relay buffer size B and packet generating 
rate A (A < p), the expected delivery delay E{iA} and the 
expected E2E delay E{T} of a packet are determined as 


E{D} 

E{r} 


1 + (n - 2 + 4 '„^_b^a)(1 - Po) 

1 - A (n - 2 + - Po) 

pLs- \ pLs 


(16) 

(17) 


where 


,S,A = 


c* • Ps 


Proof: We focus on a packet y which is the HoL packet 
of the source queue at time slot t, then in the next time slot, 
y will be delivered to its destination with probability psd, be 
forwarded to a relay node with probability psr ■ (1 — Po), and 
still stay in the source queue with probability 1 —Ms- Thus, the 
delivery process of packet y can be modeled as an absorbing 
Markov chain as illustrated in Eig.|5] where S, R and D denote 
the states that y is in source queue, forwarded to a relay, and 
delivered to its destination, respectively. We denote by Xs 
and Xu the average transition times from the transient states 
S and R to the absorbing state D, respectively. Then we have 


Xs = 1 + Xs ■ (1 — Ms) + Xr ■ Psr(l — Po), (18) 

Xs = ^+^R-Psr{l-Po) ^ 

Ms 

We denote by a probability vector P = (pojPij • • • ,Pb-i) 
the steady state distribution of the cotTesponding relay queue 
of 2 / in a relay node R, where each element pi denotes the 
probability that when y enters the relay queue, there are i 
packets already in this queue. Notice that the location of 
each node is stationary and ergodic with stationary distribution 



uniform on the network area, thus when R conducts the R-D 
transmission with probability prd in a time slot, it will deliver a 
packet for each of the n—2 traffic flows with equal probability. 
Thus, if there are i packets already in the relay queue of y, the 
expected time elapsed for y to be delivered to its destination 


is (i + 1) • 



. Then we have 


— n — 2 n — 2 n — 2 

Xr=Po ■ -l-2pi --1- \-B ■ pB-i ■ 


Prd 


Prd 


Prd 


n 


-2 


{1 + Pi + 2p2 H- \-{B - 1)pb-i} 


Prd 

n — 2 —* 

' ( 13 ” ^relay ) 5 

Prd 


( 20 ) 

( 21 ) 

( 22 ) 


where (EB follows from the normalization property of a 
probability vector, and is the average queue length of 

a relay queue, under the condition that the relay buffer is not 
full. 

We denote by H* = ,7r2^_J the limiting 

distribution of the level of a relay buffer, under the condition 
that the relay buffer is not full, then we have 


TT 


* _ 

Li — 


1 '^Lb 


C^Pl 


E B — 1 


Cfc • p 


k ’ 
s 


(23) 


and the corresponding conditional average number of packets 
occupying the relay buffer E{L*} is given by 


B-l 


E{i*}=E 


I ■ TTr. = 


2=0 




= 


,B,\- 


(24) 


Since these buffered packets are destined to each of the n — 2 
destinations with equal probability, then we have 


B 


relay 


E{L*} 

n — 2 


(25) 


Substituting the results of (1221) . (l24l i and (l25T l into ( fT^ . 
the average transition times from the transient state S to the 
absorbing state D is determined as 


Xs 


1 + (n — 2 + 'I'n,S,A)(l ~ Po) 
Ms 


(26) 


Notice that E{L>} = Xs, the result (fT^ follows, and then the 
result (fTTl) follows from E{T} = EjQ} -|- E{iA}. 


Remark 2: Similar to the two-hop scenario, in the multi¬ 
hop MANETs the packet delivery process and occupancy 
behaviors of a relay buffer can be also modeled as an AMC 
and a QBD process respectively, so it is expected that our 
proposed theoretical framework for E2E delay modeling of 
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Fig. 6. Illustration of the application of our theoretical framework. 


two-hop MANETs can be also helpful for that of the multi¬ 
hop scenarios. It is notable, however, the state transition matrix 
of the QBD process will be different under the two scenarios, 
and also that with the multi-hop network scenarios there will 
be multiple transient states in the AMC. 

Based on Theorem[Tl we can further extend our delay results 
to the buffer-unlimited scenario (i.e., B —?> oo), which is shown 
in the following corollary. 

Corollary 1: Considering the relay buffer size tends to 
infinity (B oo), then E{Z3} and E{r} are determined as 

nm = — i (27) 

B^oo Psd Psr Psd Psr ^ 

nn = . ( 28 ) 

B—>oo Psd + Psr ^ 

Proof: See Appendix |0 for the proof. ■ 

Remark 3: Notice that when B ^ oo, the results of 
Lemma 12] and Corollary [T] are coincident with the capacity and 
delay derived in [jS], where the relay buffer size is assumed to 
be infinite. 


V. Case Studies 

In this section, we conduct case studies to illustrate the 
application of our theoretical framework for the E2E delay 
modeling in buffer-limited MANETs. As illustrated in Eig. |6| 
for a given network scenario, the corresponding psd, Psr 
and Prd should be determined first, then with the inputs of 
these probabilities, by sequentially executing the relay buffer 
analysis module and delay analysis module, this framework 
finally returns the delay results. The details of the application 
of our framework under specific network scenarios are shown 
in the following subsections. 

A. Cell-partitioned MANET with LS-MAC 

We first consider a cell-partitioned MANET with local 
scheduling based MAC protocol (LS-MAC) ||8]. The whole 
network area is partitioned into m x m non-overlapping cells 
of equal size. In a time slot, each cell can support only one pair 
of nodes for packet transmission, concurrent transmissions in 
different cells will not interference with each other, and nodes 
within different cells cannot communicate. At the beginning 
of each time slot, all nodes in a cell contends for the wireless 
channel access using a DCE-style mechanism ll35l . In addition 
to these network settings, the MANET also meets the set of 
assumptions described in Section III-AI 

With the detailed information of network settings, we then 
determine the corresponding probabilities psd, Psr and prd, 
provided in the following lemma. 




(a) Transmission range of a node. (b) Illustration of EC-MAC schedul¬ 
ing. 

Fig. 7. A cell-partitioned MANET with power control and EC-MAC 


Lemma 3: For the concerned cell-partitioned MANET with 
LS-MAC, the probabilities psd, Psr and prd are given by 


m 

Psd = - 

n 


m 


n — 1 


— 1 ,m 
+( 


-1 




- 1 


1 


77—1 




(29) 


Psr —Prd 

1 I" rri^ — 1 
“2 \ n-1 





n—1 


(30) 


Proof: See Appendix iDl for the proof. ■ 

Given the number of nodes n and relay buffer size B, 
substituting formulas (|29] | and (l30l) into formula (fTTl l. we 
first determine the throughput capacity p of such a MANET. 
Then with any packet generating rate A < /i, we substitute 
formulas (|2^ and (l30l) into equation dSil to determine the 
corresponding ROP Po, and ps can be further determined by 
formula ([T]i. Substituting A, Po and ps into formulas (fTSt . 
(fThl l and (fTTl i. we finally obtain the results of queuing delay, 
delivery delay and E2E delay respectively for the concerned 
buffer-limited MANET. 


B. Cell-partitioned MANET with Power Control and EC-MAC 
We then consider a more general cell-partitioned MANET 
which applies the power control and Equivalent-Class based 
MAC protocol (EC-MAC) IB], HSj, As shown in 
Eig. |7(a)| the transmission range of a transmitter TX covers 
a set of cells which have a horizontal and vertical distance of 
no more than v —1 cells away from its own cell. Meanwhile, 
as illustrated in Eig. |7(b)| all cells are divided into different 
ECs, where any two cells in the same EC have a horizontal 
and vertical distance of some multiple of e cells. Thus, the 
MANET contains in total ECs and ECs are activated 
alternatively as time evolves. Suppose that at time slot t, a 
transmitter TXq in an active cell will transmit a packet to its 
receiver RXq, in order to ensure the transmission successful, 
according to the Protocol Model ll36l it should satisfy that 

dTXi,RXo > (1 + ^)dTXo,RXo, ( 31 ) 

where TXi denotes a concurrent transmitter in any one of the 
other active cells, dij denotes the distance between nodes i 
and j, and A is a guard factor. Thus we have 

e — o > (1 -\- A)-\/2i/, 


(32) 























































































and e should be set as 


1 


-Theoretical 

O I.i.d,simulation 
^ Random walk, simulation 


£ = min{|’(l + A)-\/2:/+ (33) 


Regarding the corresponding probabilities Psdi Psr and pr-d 
of this type of MANETs, we have the following lemma. 

Lemma 4: For the concerned MANET with power control 
and EC-MAC, the probabilities psd, Psr and prd are given by 



r-^ 


—1 —(E—l)n 1 


n — 1 


i{n — 1) 




(34) 


Psr —Prd 

1 


m? — T 


n—1 



(35) 


where E = {2iy — 1)^. 

Proof: See Appendix [D] for the proof. ■ 

By applying the same operations of our theoretical frame¬ 
work as the previous subsection, we can obtain the delay 
results for the concerned MANETs. 


C. Other Network Scenarios 

Notice that to apply our theoretical framework for the E2E 
delay modeling, it only needs to determine the inputs of the 
framework, i.e., the probabilities psd, Psr and prd- Thus, this 
framework also has the great potential to be applied to many 
other network scenarios. For example, for the MANETs where 
the 2HR routing scheme is administered by cell il (not by 
each node in our case), and the ratio between S-R and R-D 
transmission can be changed ifTSlI (we fix the ratio as 0.5), 
these probabilities can be determined. Recently, Chen et al. 
iTT^ has reported that how to compute these probabilities for 
a MANET under the continuous network model and the Aloha 
MAC protocol, then with these probabilities as the inputs of 
our framework, the corresponding delay analysis under the 
practical buffer constraint can be conducted. 



0.3 0.4 0.5 0.6 

Workload ^jj.) 


(a) ROP under LS-MAC. 



Fig. 8. Theoretical and simulated ROP performance. Case 1: n = 32, m = 

4, S = 5. Case 2: n = b0,m = 6, B = b. 


B. Validation 


VI. Simulation Results 

In this section, we hrst conduct simulations to validate our 
theoretical framework for the E2E delay modeling in buffer- 
limited MANETs, then provide discussions about the impacts 
of network parameters on delay performance. 

A. Simulation Settings 

For the validation of our theoretical framework and delay 
results, a specihc C-n- simulator was developed to simulate 
the packet generating, queuing and delivery processes in a 
cell-partitioned MANET llJTl . where the network settings, 
including relay buffer size B, number of nodes n, partition 
parameter m, packet generating rate A and the mobility model 
can be flexibly adjusted to simulate the network performance 
under various scenarios. For the network scenario with power 
control and EC-MAC, we set = 1 and A = 1 Il38l . The 
duration of each task of simulation is set to be 2 x 10® time 
slots, and we only collect data from the last 80% of the time 
slots in each task (the system will be in the steady state with 
high probability), to ensure the accuracy of simulated results. 


First, we provide plots of the theoretical and simulated 
ROP performance under two network scenarios in Fig. [8] 
and for each scenario we consider two cases (case 1; n = 
32, m = 4, i? = 5, and case 2: n = 50, m = 5, B — 5) 
and two mobility models (the i.i.d mobility model and the 
random walk model). The workload is dehned as A//i. We 
can see from Fig. 0 that the simulation results match nicely 
with the theoretical ones for all the cases, which indicates that 
our framework is highly efficient in depicting the occupancy 
behaviors of the relay buffer in buffer-limited MANETs. 

Then, with the same network settings, we provide plots of 
the theoretical and simulated E2E delay results in Fig. |9] 
It is observed from Fig. |9] that all the simulation results 
can match the corresponding theoretical curves very nicely, 
indicating that; 1) our theoretical framework is highly efficient 
for the E2E delay modeling in buffer-limited MANETs; 2) the 
framework is very general since it can be applied to various 
network scenarios. Another observation of Fig. |9] is that the 
packet E2E delay increases sharply as the packet generating 
rate A approaches a specihc value (e.g., under LS-MAC and 
case 1, the value is around 0.038), which serves as an intuitive 










9 



(a) End-to-end delay under LS-MAC. 


O Li.d, simulation 
* Random walk, simulation 


LS-MAC 


-o- 


* 

^ o 

* 

o 


\ 


o 

* * ^ ^ 
0 . 0 ^ 


Power control and EC—MAC 


0.4 0.6 0.8 1 

Workload ^|i) 



(b) End-to-end delay under power control and EC-MAC. 


Fig. 9. Theoretical and simulated end-to-end delay performance. Case 1: 

n = 32, m = A, B = 5. Case 2: n = 50, m = b, B = 


Fig. 10. Per node achievable throughput. LS-MAC: n = 200, m = 10, B = 
5. Power control and EC-MAC: n = 32, m = A, B = 10. 



Fig. 11. Delivery delay vs. workload (X/fi) under different settings of relay 
buffer size, n = 32, m = 4. 


impression of its corresponding throughput capacity fi. 

To further validate our framework on throughput capacity. 
Fig. [TO] summarizes the simulation results on the achiev¬ 
able per node throughput, where two network scenarios with 
different throughput capacity (LS-MAC; n — 200,m = 
10, B = 5,/rLS = 6.5 x 10“^) and (power control and EC- 
MAC; n = 32, m = A, B = 10, fj,EC = 3.3 x 10“^) are 
presented. We can see that for each network scenario there, 
the per node throughput first monotonously increases before 
the workload reaches 1, and then remains a constant which is 
just the corresponding throughput capacity when the workload 
exceeds 1. Thus, it indicates that the proposed framework is 
also efficient in depicting the per node throughput capacity 
behavior of the buffer-limited MANET. 

C. Performance Discussion 

With the help of our theoretical framework for the E2E 
delay modeling, we explore how the network parameters 
affect the the delay performance of a buffer-limited MANET. 
Without loss of generality, we consider here a cell-partitioned 
MANET with LS-MAC. 


We first summarize in Fig.[TT]that how the expected delivery 
delay E{Z?} varies with the workload. A very interesting 
observation is that under the buffer-limited scenarios (B = 5 
and B = 20), as workload increases, E{L)} first increases to 
a maximum and then decreases. This is due to the reason that 
the effects of workload on E{L)} are two folds. On one hand, 
a larger workload will lead to a longer relay queue length, 
which further leads to a higher delay in the relay queue; on 
the other hand, a larger workload will lead to a higher ROP, 
which further leads to a lower probability that a packet to be 
delivered by a two-hop way, such that E{L)} decreases. Since 
the latter effect, the delivery delay under a small relay buffer 
is lower than that under a large one. 

Fig. [TOjshows the relationship between the expected E2E de¬ 
lay E{T} and packet generating rate A. We can see that under 
buffer-limited scenarios, as A increases, E{T} doesn’t increase 
all the time because the delivery delay will decrease when 
A exceeds a specific value, however when A approaches the 
corresponding throughput capacity, E{T} increases sharply 
because the queuing delay tends to infinity. It also can be 
seen that when A is small, E{T} under B = 5 is smaller than 
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Fig. 12. End-to-end delay vs. packet generating rate A under different settings 
of relay buffer size, n = 32, m = 4. 



Fig. 13. End-to-end delay vs. relay buffer size, n = 32, m = 4. 

that under B = 20, since both of the queuing delay under 
two settings are small, but a small relay buffer can lead to a 
small delivery delay. However, with A getting larger and larger, 
E{T} under B = 5 hnally exceeds that under B = 20, and 
tends to inhnity earlier. It indicates that increasing the relay 
buffer size can ensure the E2E delay limited for a larger region 
of packet generating rate. 

We illustrate in Eig. [13] how E{T} varies B under the 
settings of (n = 32, m = 4, A = {0.01,0.02}). According 
to formula (ITSl i. p. = 0.0227 when B = 1, and n increases as 
B increases. Thus, for A = 0.01 which is much smaller than 
0.0227, E{T} increases as B increases and hnally tends to 
a constant 206.92 which can be determined by formula (l28T l. 
While for A = 0.02 which is very close to the ^ under B = 1, 
E{T} under B = 1 is very large. With B increasing, the 
corresponding /r increases, leading to the E{T} hrst decreases, 
then increases and hnally tends to a constant 221.65. 

We further illustrate in two 3D hgures (Eig. [T4| and Eig. [TsT i 
that how E{r} is inhuenced by (n. A} and {n,B}, respec¬ 
tively (the ratio of n to the number of cells keeps as 2). We 
can see that the variations of E{r} with n are complicated, 
but in general E{T} increases as n increases. A more careful 
observation is that when n increases, E{T} hrst increases 



Fig. 14. End-to-end delay vs. packet generating rate A and number of nodes 
n. S = 5. 



Fig. 15. End-to-end delay vs. relay buffer size B and number of nodes n. 
A = 0.02. 

almost linearly when A is much smaller than fji, then increases 
quickly when A approaches /r. Eor example, these behaviors 
can be found in Eig. fT4l under A = 0.23 and in Eig. [TSl under 
B = l. 

VII. Conclusion 

This paper represents a signihcant step towards the ex¬ 
act end-to-end delay modeling of practical buffer-limited 
MANETs. With the help of the theories of Eixed-Point, QBD 
process, EMC and AMC, a novel theoretical framework has 
been developed to efficiently depict the highly dynamics in 
such networks. This framework is general in the sense that 
it can be applied to the MANETs with any mobility model 
that leads to the uniform distribution of the locations of 
nodes, and any MAC protocol as long as the probabilities 
Psd, Psr and Prd there can be determined. Also, it is expected 
the framework can shed light on the E2E delay modeling 
for multi-hop MANETs. Extensive simulations have been 
conducted to validate the efficiency and applicability of our 
framework, and some interesting theoretical hndings about 
the impacts of network parameters on delay performance have 
been discussed. 

Notice that our E2E delay modeling for buffer-limited 
MANETs is based on the 2HR routing and relay buffer 
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constraint, so one of the future research directions is to extend 
our study to the E2E delay modeling for MANETs with 
multi-hop routing schemes and more practical buffer constraint 
on both source buffer and relay buffer. Another appealing 
direction is to explore the delay modeling for the concerned 
MANETs with the consideration of more practical network 
settings (such as the wireless channel fading) and apply the 
well-known NS-2 network simulator for model validation. 


Appendix A 
Proof of Lemma[T] 

We denote by A^j the packet aiTival rate at the relay buffer in 
S. Due to the ergodic and uniform properties of node mobility, 
for node S, each of the remaining n — 2 nodes (except D) is 
likely to serve as its relay with equal probability. Due to the 
symmetry of nodes and traffic flows, for S serving as a relay, 
all the other n — 2 nodes are likely to forward packets to its 
relay buffer. Notice that the packet departure process of source 
queue is a Bernoulli process with rate A, and from ([T]i we can 
see the ratio of S-R transmission to the whole packet departure 

is jhen we have 

ij-bW 

\ f \ Psr(l~Po(A)) ii ^ 

\R = [n- 2)A--/(n - 2) 

= Ps{X)psr (1 - PoW) ■ (36) 

Regarding the transition probability from state (i,j) to the 
states in its adjacent upper level L^+i, we have the following 
equation 


P(^,j),Li+i • (1 “ Po{X)) + 0 • po(A) = Afi. (37) 
Combining dJTl l with (l36l l we have 


Substituting 


P{i,j),Li^i — Ps(A)ps 

into © we have 

^ p^(^X)psr 


(38) 


Regarding the transition probability from state {i,j) to the 
states in its adjacent lower level Li_i, when S conducts a R- 
D transmission with probability prd, due to the ergodic and 
uniform features of node mobility, it will choose one of the 
other n — 2 nodes as its receiver with equal probability. Then 
we have 


P{^,3)M-l=Prd■ (39) 

To determine the conditional probability Pj\Li, we utilize 
the following Occupancy technique 1391 . Considering the relay 
buffer on level i, where each of these i buffered packets may 
be destined for any one of the other n — 2 nodes (except S 
and D), the number of all possible cases is determined 
as 


Nl, 




(40) 


Considering the condition that these i packets are destined 
for only j different nodes, then the number of possible cases 
A(j is determined as 


^(ij) ~ 


n — 2 

j 


(j -!) + (*- j) 
i-j 


Due to the ergodic and uniform features of node mobility, each 
of these cases occurs with equal probability. According to the 
Classical Probability, Pj\Li is then determined as 


Pj\Li - 


N, 


Nl. ("-f+0 


(ij) 


(42) 


It can be easily verified that ^ Pj\Li = 1- Combining ( l39b . 
(l42l i and © we have 


(M-i) A f ("J7 ■ (}:1) i 

Pi =E|—ps+Jj—^ f''!— 

_ Prd ^r/n-3A /i-lAI 

” (””^0 ’ ^ H ^ / ' V J Jj 

/'n— 4 +i\ 

(#7° (43) 


Appendix B 
Proof of Lemma|2] 

Since the relay buffer is limited, it is always stable. Thus, we 
only need to consider the stability of source queue. The source 
queue is a Bernoulli/Bernoulli queue with packet generating 
rate A and a corresponding service rate Ps{X). Based on the 
queuing theory, the source queue is stable (resp. unstable) 
when A < Ps{X) (resp. A > /is(A)). 

As A increases, Po(A) increases since the network is more 
congested, and /is (A) decreases according to formula ©. 
Thus, from the definition of throughput capacity, /i should 
satisfy that 

/i = A* = /is(A*). (44) 


It is notable that when A approaches A*, lim Ps{X) = 1, and 
from © we have 


PoiX*) = 


Cb 


J 2 k =0 


n — 2 
n — 2 + B 


(45) 


Then p is determined by substituting (l45l l into ©. 


Appendix C 

Proof of Corollary[T] 


We denote by F{ps), G{ps) the sums of infinite series 
Si>0 Cip* and Y.i>Q iCipl, respectively. Notice that F{ps) 
is the Taylor series expansion from (1 — then we have 


FiPs) = 


(1-Ps) 


n-2 ’ 


G{ps) = Ps ■ F'{pa) = {n-2) 




Eurther we have 


and 


,A = 


GiPs) , Ps 


F{Ps) 


= {n-2) 


I - Ps 


Po = lim Cb • pf • (1 - Ps)" ^ 

B —¥00 B —¥00 

< lim {B + n)'^pf < lim 2"B"pf 

B —¥00 B —¥00 

Tl^ 

= lim ^= 0, 


B^oo (— In Ps)' 


(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


(41) 
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where ( ISTT l is obtained by utilizing the L’Hopital’s rule recur¬ 
sively. 

Substituting (l48l l and ( fSTT l into Theorem [T] we can obtain 
dZTl i and (l28T l directly. 


Appendix D 

Proofs of Lemma[3]and Lemma|4] 

For a cell-partitioned MANET with LS-MAC, the event that 
node S conducts a S-D (resp. S-R or R-D) transmission in a 
time slot can be divided into the following sub-events: (1) D 
is (resp. is not) in the same cell with S; (2) other k out of 
n — 2 nodes are in the same cell with S, while the remaining 
n — 2 — k nodes are not in this cell; (3) S contends for the 
wireless channel access successfully. Thus we have 



For a cell-partitioned MANET with power control and 
EC-MAC, by applying the similar approach and algebraic 
operations we have 



and 

Psr — Prd 


1 m? — r 

2 ^^ 
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